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Abstract. In this paper, we develop an exact test for testing 
exponentiality against renewal increasing mean residual life class. 
Pitman’s asymptotic efficacy value shows that our test perform 
well. Some numerical results are presented to demonstrate the 
performance of the testing method. We also discuss how the pro¬ 
posed method incorporates the right censored observations. 
Keywords: Exponential distribution; Renewal increasing mean resid¬ 
ual life; Replacement model; Shock Model; U-statistics. 


1. Introduction 

When a device is experiencing random number of shocks governed by a 
homogeneous Poisson process, the concept of renewal increasing mean 
residual life is very much useful to study the properties of age replace¬ 
ment model. In this context, test for exponentiality against renewal 
increasing mean residual life class is used to determine whether to 
adopt a planned replace model over unscheduled one. Sepehrifar et 
al. (2015) developed a non-parametric test against RIMRL s /i 0cfc class 
and obtained a critical region based on the asymptotic theory of U- 
statistics. We noted that the critical region developed by Sepehrifar et 
al. (2015) is incorrect (see Remark 2.1). Motivated by Sepehrifar et 
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al. (2015), we develop an exact test for testing exponentiality against 
RIMR L s /, oc /. class. We also obtain the correct critical region of the 
asymptotic test proposed by Sepehrifar et al. (2015). The case with 
censored observations also addressed. 

The rest of the paper is organized as follows. In Section 2, we pro¬ 
pose an exact test for testing exponentiality against RIMRLs/joc*. class 
and then calculate the critical values for different sample sizes. The 
asymptotic normality proposed test statistic is proved in Section 3. 
The Pitman’s asymptotic efficacy value is also given in this section. In 
Section 4, we report the result of simulation study carried out to assess 
the performance of the proposed test. In Section 5, we discuss how to 
incorporate right censored observation in our study. 

2. Exact Test 

Let X be the lifetime of a device which has absolutely continuous distri¬ 
bution function F(.). Suppose F(x) = P(X > x) denotes the survival 
function of X at x. Also let /i = E(X) = J 0 °° F(t)dt < oo. Assume 
that the device under consideration is experiencing a random shock. 
Suppose N(t) denotes the total number of shocks up to time t with 
probability mass function P(N(t) = j ) = F j (t) — F j+1 (t),j = 0,1, 2,.... 
Suppose that the random variable Wj,j = 0,1, 2, quantify the amount 
of hidden lifetime absorbed by the jth shock with Wq = 0 and hav¬ 
ing common distribution function G(x) = P{Wj < x). The total 
cumulative life damage up to time t is defined as Z(t) = Wj 

with the cumulative distribution function Q(x) = P(Z(t) < x) — 
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£~ 0 G^^F^(t) — F^ +1 \t)]. It is assumed that the unit fails when 
the total life-damage exceeds a pre-specihed level x > 0. We refer to 
Glynn and Whitt (1993), Roginsky (1994) and Sepehrifar et al. (2015) 
for discussion related this framework. 


Let X* = X — Z(t) be the residual lifetime of an operating device 


with cumulative damage Z{t). Note that the realizations of X* is 
available to us for further analysis. Consider a device subjected to 
N(t) number of shocks up to time t. Given that such a device is in 
an operating situation at time instant t after installation, the MRL 
function of X* denoted by m*{t ) is defined by m*(t) = E(X* — t\X* > 
t ). Note that the total life-damage will not exceed the threshold level 
x. From the definitions it is evident that the random variables X* and 
Z(t ) are independent. Next we give the definition of RIMRL s /i OC fc class 
(Sepehrifar et ah, 2015). 

Definition 2.1. The mean residual life of a device under shock model 
(MRL shock ) at time t is defined as 



where r(z) = f* F(z + w)dQ(w). 

Definition 2.2. The random variable X belongs to the RIMRL s ho C k 
class if the function m*(t ) is a non-decreasing function for all t > 0. 

We are interested to test the null hypothesis 


H o : F* is exponential 
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against the alternatives 

Hi : F* is RIMRL s / loc fc (and not exponential), 


on the basis of a random sample X*, X %,..., X*] from an absolutely 
continuous distribution function F*. For the above testing problem 
Sepehrifar et al. (2015) proposed a non-parametric test based on the 
departure measure A*(F*) defined by 


A *(F*) 




A (F*) 


where A(F*) = Ef*(min(X(, X 2 ) — ^X() and p* = E(X)). Based on 
U-statistics theory Sepehrifar et al. (2015) obtained the following test 
statistic 


A* = 


A 

X*’ 


( 1 ) 


where X* = ± EHi x t and ^ = F(hi) ELi EJ<iy=i H X *,X*) with 
h(Xl, X* 2 ) = min(X*, X 2 ) — |A,*. Hence the test procedure is to reject 
the null hypothesis H 0 in favour of Hi for large values of A*. 


Remark 2.1. For the testing problem discussed here, Sepehrifar et al. 
(2015) obtained a critical region based on the asymptotic variance ■J-. 
However as we shown (see Section 3) the asymptotic variance is X. We 
could not find enough details in their paper to explain the discrepancy. 

Motivated by this discrepancy, next we develop an exact test based 
on the test statistics A* and calculate the critical values for different 
sample size. We use a result due to Box (1954) to find the exact null 
distribution of the test statistic. 
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Theorem 2.1. Let X* be continuous non-negative random variable 
with F*(x) = e - i. Let X*. Xtf ,.... X* be independent and identical 
samples from F*. Then for fixed n 

n n , 

P(A* > X) = ^ II (' rf l,U _f ) I ( X ’ d i,n )» 

. ' tti.n U'j.n ' 

1=1 

provided d i)U ^ dj >n for i ^ j, where 


I{x,y) 


1 if x < y 
0 if x > y 


and 


(n — 2i + 1) 
“ 2 (n - 1) 


Proof: First we express the test statistics in terms of order statistics. 
Note that 


0 n n ~ ,v 


nin — 1) , , . . , 

z=i ]<%\]=i 


n\ 


i{n - 1) 

where XT , i = 1,2, is the z-th order statistics based on the ran¬ 
dom sample -Xjf, X|, ...X*; from F*. After some algebraic manipula¬ 
tion, we can express the equation (pQ) as 

M^ryEILi (3n-4z+l)X ( * } 


A* = 


X* 


( 2 ) 


Rewrite the denominator of the equation ([2]) as 


* = £*<*> 


~(n — i + l) 2 
. n(n — 1) 



n(n — 1) 


(n + 1) ' 

2 n(n — 1). 


Or 


A 


n _ \ ( n ~ i + 1 ) 2 

n - 1) ^ (i) i n 2 

' i =1 


(n — z) 2 (n + 1)' 
n 2 2 n 2 . 
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Hence, in terms of the normalized spacings, Di = (n — i + l)(X?h — 
^(*_i))’ with X* = 0, we can express the test statistics as 


A* 


Ek d i,nDj 

E”= 1 A 


where d^s are given by 


di 


(E T) [( "- i + 1) 

(n — 2z + 1) 

2(n — 1) 


A + l) i 

2 J 


Note that the exponential random variable with rate | is distritbuted 
same as the x 2 random variable with 2 degrees of freedom. Hence the 
result follows from Theorem 2.4 of Box (1954). 

The critical values of the exact test for different n are tabulated in 
Table 1. 


Table 1. Critical values of the exact test 


n 

90% level 

95% level 

97.5% level 

99% level 

2 

0.4000 

0.4500 

0.4750 

0.4900 

3 

0.2764 

0.3419 

0.3883 

0.4292 

4 

0.2189 

0.2678 

0.323 

0.3693 

5 

0.1883 

0.2383 

0.28 

0.325 

6 

0.1679 

0.2131 

0.2508 

0.2927 

7 

0.1529 

0.1944 

0.2293 

0.2682 

8 

0.1413 

0.1799 

0.2125 

0.2492 

9 

0.1319 

0.1682 

0.1989 

0.2336 

10 

0.1243 

0.1586 

0.1877 

0.2208 

15 

0.0993 

0.1271 

0.1508 

0.178 

20 

0.0852 

0.109 

0.1295 

0.1531 

25 

0.0758 

0.097 

0.1153 

0.1363 

30 

0.0689 

0.0882 

0.1049 

0.1241 

40 

0.0594 

0.0761 

0.0905 

0.1072 

50 

0.0529 

0.0679 

0.0808 

0.0957 

75 

0.0431 

0.0552 

0.0658 

0.078 

100 

0.0373 

0.0477 

0.0569 

0.0675 
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3. Asymptotic properties 

In this section, we prove asymptotic normality of the proposed test 
statistic. Making use of the asymptotic distribution we also calculate 
the Pitman’s asymptotic efficacy of the test. As mentioned Sepehrifar 
et al. (2015) showed that the test statistic has limiting normal distribu¬ 
tion, however they incorrectly stated the asymptotic variance. Hence 
we give the following results to correct the error occurred in their study. 

Theorem 3.1. The distribution of y/n( A — A(F*)), as n —>• oo, is 
Gaussian with mean zero and variance 4 af, where a\ is the asymptotic 
variance of A and is given by 

af = ±Var(2X*F*(X*) + 2^ ydF\y) - ^X*). (3) 

Corollary 3.1. Let X* be continuous non-negative random variable 
with F*(x) = e _ y then the distribution of y/nA, as n —?• oo, is Gauss¬ 
ian with mean zero and variance a^ = y^. 

Corollary 3.2. Let X* be continuous non-negative random variable 
with F*{x) = e~%, then the distribution of y/nA* , as n —» oo, is Gauss¬ 
ian with mean zero and variance a% = 

Apart from the exact test we can construct an asymptotic test based 
on the asymptotic distribution of A*. Hence in case of the asymptotic 
test, for large values of n, we reject the null hypothesis H 0 in favour of 
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the alternative hypothesis H \, if 

a/I2to(A*) > Z a , 


where Z a is the upper a-percentile of N(0, 1). In fact, this is the correct 
critical region of the test proposed by Sepehrifar et ah (2015). 

Next we study the asymptotic efficiency of the test. The Pitman’s 
asymptotic efficacy is the most frequently used index to make a quan¬ 
titative comparison of two distinct asymptotic tests for a certain sta¬ 
tistical hypothesis. The Pitman’s asymptotic efficacy (PAE) is defined 
as 

PAE( A*) = ldX A )I A ~* A ° 


where Ao is the value of A under Hq and cTq is the asymptotic variance 
of A* under H 0 . In our case, the PAE is given by 


PAE(A*(F*)) = 


<7o 


= vd2(W'(A„) - WdA„K'(A„)), 


where W = E(min(Xl, X%)) and p* is the mean of X* under the 
alternative hypothesis and the prime denotes the differentiation with 
respect to A. We calculate the PAE value for three commonly used 
alternatives which are the members of RIMRL s / )OC fc class 

(i) the Weibull family: F*(x) = e~ xX for A > 1, x > 0 

(ii) the linear failure rate family: F*(x) = e^~ x ~^ x2 ^ for A > 0, x > 0 

(iii) the Makeham family: F*{x ) = e~ x ~ x ^ e x+x ~ l ) for A > 0, x > 0. 
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By direct calculations, we observe that the PAE for Weibull distribu¬ 
tion is equal to 1.2005; while for linear failure rate distribution and the 
Makeham distribution these values are, 0.8660 and 0.2828, respectively. 

Next we compare the performance of the proposed test with some 
other tests available in the context of age replace model by evaluating 
the PAE of the respective tests. We compare our test with that tests 
proposed by Li and Xu (2008) and Kayid et al. (2013). The Table 
2 gives the PAE values for different test procedures. From the Table 
2, it is clear that our test is quite efficient for the Weibull and linear 
failure rate alternatives. Note that the test proposed by Kayid et al. 
(2013) has good efficacy for Makeham alternative even though their 
test shows poor performance against the other two given alternatives. 

Table 2. Pitman’s asymptotic efficacy (PAE) 


Distribution 

Proposed test 

Li and Xu (2008) 

Kayid et al. (2013) 

Weibull 

1.2005 

1.1215 

0.4822 

Linear failure rate 

0.8660 

0.5032 

0.4564 

Makeham 

0.2828 

0.2414 

2.084 


4. Simulation study 

Here we report a simulation study for evaluating the performance of 
our asymptotic test against various alternatives. The simulation was 
done using R program. 

First we fold the empirical type 1 error of the proposed test. Since 
the test is scale invariant, we simulate random sample from standard 
exponential distribution. The simulation is repeated for ten thousand 
times with different values of n and is reported in Table 3. From 
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Table 3. Empirical type 1 error of the test 


n 

Type 1 Error (5% level) 

Type 1 Error (1% level) 

10 

0.0635 

0.0123 

20 

0.0540 

0.0115 

30 

0.0518 

0.0107 

40 

0.0520 

0.0110 

50 

0.0517 

0.0107 

60 

0.0516 

0.0102 

70 

0.0515 

0.0102 

80 

0.0511 

0.0100 

90 

0.0504 

0.0103 

100 

0.0504 

0.0104 


the Table 3 it evident that the empirical type 1 error is a very good 
estimator of the size of the test even for small sample size. 

For Ending empirical power against various alternatives, we simulate 
observations from Weibull, linear failure rate and Makeham distribu¬ 
tions with different values of A where the distribution functions were 
given in the Section 3. The empirical powers for the above mentioned 
alternatives are given in Tables 4, 5 and 6. From these tables we can 
see that empirical powers of the test approaches one when the 9 values 
are going away from the null hypothesis value as well as when n takes 
large values. 

TABLE 4. Empirical Power: Weibull distribution 



A = 1.2 

A = 1.4 

A = 

1.6 

A = 1.8 

n 

3Vo 

1% 

37o 

1% 

5% 

1% 

3Vo 

1% 

60 

0.50 

0.23 

0.93 

0.76 

0.99 

0.97 

1.00 

0.99 

70 

0.55 

0.27 

0.96 

0.84 

0.99 

0.99 

1.00 

1.00 

80 

0.60 

0.31 

0.98 

0.89 

0.99 

0.99 

1.00 

1.00 

90 

0.64 

0.36 

0.99 

0.93 

0.99 

0.99 

1.00 

1.00 

100 

0.69 

0.41 

0.99 

0.95 

1.00 

0.99 

1.00 

1.00 
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Table 5. Empirical Power: Linear failure rate distribution 



A = 0.2 

A = 0.4 

A = 

0.6 

A = 0.8 

n 

5% 

1% 

5% 

1% 

5% 

1% 

5% 

1% 

60 

0.49 

0.22 

0.68 

0.38 

0.79 

0.51 

0.87 

0.65 

70 

0.55 

0.27 

0.74 

0.46 

0.84 

0.61 

0.92 

0.74 

80 

0.60 

0.32 

0.80 

0.53 

0.89 

0.68 

0.94 

0.81 

90 

0.65 

0.36 

0.83 

0.59 

0.91 

0.74 

0.97 

0.86 

100 

0.69 

0.41 

0.87 

0.65 

0.94 

0.80 

0.98 

0.90 


Table 6. Empirical Power: Makeham distribution 



A = 0.2 

A = 0.4 

A = 0.6 

A = 0.8 

n 

5% 

1% 

5to 

1% 

5% 

1% 

5% 

1% 

60 

0.37 

0.14 

0.49 

0.22 

0.65 

0.36 

0.87 

0.63 

70 

0.42 

0.17 

0.55 

0.26 

0.72 

0.43 

0.92 

0.72 

80 

0.46 

0.20 

0.60 

0.31 

0.78 

0.49 

0.94 

0.79 

90 

0.51 

0.23 

0.65 

0.35 

0.82 

0.56 

0.96 

0.84 

100 

0.55 

0.27 

0.70 

0.40 

0.86 

0.62 

0.98 

0.90 


5. The case of censored observations 


Next we discuss how the censored observations can be incorporated 
in the proposed method. Suppose we have randomly right-censored 
observations such that the censoring times are independent of the life¬ 
times. Under this set up the observed data are n independent and 
identical copies of (Y*, S), with Y* = min(X *, (7), where C is the cen¬ 
soring time and 5 = I(X* < C). Now we need to address the testing 
problem mentioned in Section 2 based on n independent and identical 
observation {(Y^*, Si), 1 < i < n}. Observe that S t — 1 means i th object 
is not censored, whereas Si = 0 means that i th object is censored by 
C , on the right. Usually we need to redefine the measure A*(F*) to 
incorporates the censored observations. The U-statistics formulations 
helps us to solve the problem in an easy way. Using the right-censored 
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version of a U-statistic introduced by Datta et al. (2010) an estimator 
A (F*) with censored observation is given by 




provided K C (Y ■*), K C (Y*) > 0, with probability 1 and K c is the Kaplan- 
Meier estimator of K c , the survival function of the censoring variable 
C. Similarly an estimator of fi* is given by (Zhao and Tsiatis, 2000) 





Hence in right censoring situation, the test statistics is given by 



and the test procedure is to reject H 0 in favour of H i for large values 


Next we obtain the limiting distribution of the test statistic. Let 
N°(t) = I(Y* <t,Si — 0) be the counting process corresponds to the 
censoring variable for the i th individual, Zj(t) = I(Y* > t). Also let A c 
be the hazard rate of C. The martingale associated with this counting 
process is given by 
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Let G(x,y) = P(X* < x,Y* < y, S = 1), x G X, H(t ) = P{Y* < t) 
and 


w{t) = 


hAx) 


m 


;I(y > t)dG(x, y), 


IXx [0,oo) K C {y~) 

where hi(x) = Eh(x, X^)- Next result follows from Datta et al. (2010) 
for the choice of the kernel h(Y*, Yf) = |(4 Y*I(Y* < Yf) + 4YfI(Y2 < 


*?) - y* - y;). 


Theorem 5.1. If Eh 2 (Y*,Y*) < oo, f Xx[0 oo) §^dG(x, y) < oo and 
J 0 °° w 2 (t)\ c (t)dt < oo, then the distribution of ■ v /n(A c — A(F*)), as 
n —$■ oo, is Gaussian with mean zero and variance 4 a\ c , where a\ c is 
given by 

*1 = Var ( h F^(Y’-) + / WMWW). 

Corollary 5.1. Under the assumptions of Theorem, \5.1\ if E(Yf) < oo, 
the distribution of y/n(A* c — A*(F*)), asn —>■ oo, is Gaussian with mean 
zero and variance 4a 2 , where 


a 


2 

c 


a 


T 


lc 

*2 


(4) 


Corollary 5.2. Let X* be continuous non-negative random variable 
with F*(x) = e _ W Under the assumptions of Theorem \5.1[ if E[Yf) < 
oo, the distribution of ^/nA*, as n —>• oo, is Gaussian with mean zero 
and variance a 2 0 , where 


°do = —Car 


A 2 


(4F*(A*) — A*)hi 

4/l c (1?—) 


w{t)dM[{t ) 


(5) 


Hence by Corollary 5.2, we know that the y/nA* has asymptotically 
normal with mean zero and a variance that can be estimated by (5) 
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and we denote it as a c o. Hence we reject the null hypothesis in favour 
of H x , if 



CtcO 


Next we study the efficiency loss due to censoring by computing the 
efficiency of our test based on A* for uncensored model and the effi¬ 
ciency of the test based on A* for censored model. As both these tests 
have same asymptotic mean, the Pitman asymptotic relative efficiency 
(ARE) of the test based on A* with respect to the test based on A* is 
given by 

2 

e = ARE(A* C ,A*) = 

a c0 

The quantity (1 — e) can be taken as a measure of the efficiency loss 
(Lim and Park, 1993) due to censoring. From the above expression it is 
clear that the ARE value is independent of the distributions belonging 
to the family of alternative hypothesis, but depends on the distribution 
of C. Next we calculate the ARE value when the censoring variable C 


TABLE 7. Asymptotic relative efficiency 


A 

0.5 

0.4 

0.3 

0.2 

0.1 

0.05 

0.01 

ARE 

0.397 

0.433 

0.480 

0.547 

0.643 

0.700 

0.741 


has logistic distribution with distribution function F(x) = The 

ARE value for different values of A is given in Table 7. Table 7 shows 
that as A decreases, the value of ARE increases and the efficiency loss 
decreases as the value of A (the amount of censoring) becomes small. 
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